A global multilayer cloud identification with POLDER/PARASOL

The detection of multilayer cloud situations is important for satellite retrieval algorithms and for many climate related applications. In this paper, we describe an algorithm based on the exploitation of the POLarization and Directionality of the Earth’s Reflectance (POLDER) observations to identify monolayered and multilayered cloudy situations along with a confidence index. Our reference comes from the synergy of the active instruments of the A-Train satellite constellation. The algorithm is based upon a decision tree that uses a metric from information theory and a series of tests on POLDER Level-2 products. We obtain a multilayer flag as the final result of a tree classification which takes discrete values between 0 and 100. Values closest to zero (resp. a hundred) indicate a higher confidence in the monolayer (resp. multilayer) character. This indicator can be used as it is, or with a threshold level that minimizes the risk of misclassification, as a binary index to distinguish between monolayer and multilayer clouds. For almost fully covered and optically thick enough cloud scenes, the risk of misclassification ranges from 29% to 34% over the period 2006–2010 and the average confidences in the estimated monolayer and multilayer characters of the cloud scenes are 74.0% and 58.2% respectively. With the binary distinction, POLDER provides a climatology of the mono/multi-layer cloud character that exhibits some interesting features. Comparisons with the performance of the Moderate-Resolution Imaging Spectroradiometer (MODIS) multilayerflag are given

Improved information about the vertical location and extent of monolayer clouds from POLDER3 measurements in the oxygen A-band

This paper describes new advances in the exploitation of oxygen A-band measurements from POLDER3 sensor onboard PARASOL, satellite platform within the A-Train. These developments result from not only an account of the dependence of POLDER oxygen parameters to cloud optical thickness τ and to the scene's geometrical conditions but also, and more importantly, from the finer understanding of the sensitivity of these parameters to cloud vertical extent. This sensitivity is made possible thanks to the multidirectional character of POLDER measurements. In the case of monolayer clouds that represent most of cloudy conditions, new oxygen parameters are obtained and calibrated from POLDER3 data colocalized with the measurements of the two active sensors of the A-Train: CALIOP/CALIPSO and CPR/CloudSat. From a parameterization that is (μs, τ) dependent, with μs the cosine of the solar zenith angle, a cloud top oxygen pressure (CTOP) and a cloud middle oxygen pressure (CMOP) are obtained, which are estimates of actual cloud top and middle pressures (CTP and CMP). Performances of CTOP and CMOP are presented by class of clouds following the ISCCP classification. In 2008, the coefficient of the correlation between CMOP and CMP is 0.81 for cirrostratus, 0.79 for stratocumulus, 0.75 for deep convective clouds. The coefficient of the correlation between CTOP and CTP is 0.75, 0.73, and 0.79 for the same cloud types. The score obtained by CTOP, defined as the confidence in the retrieval for a particular range of inferred value and for a given error, is higher than the one of MODIS CTP estimate. Scores of CTOP are the highest for bin value of CTP superior in numbers. For liquid (ice) clouds and an error of 30 hPa (50 hPa), the score of CTOP reaches 50% (70%). From the difference between CTOP and CMOP, a first estimate of the cloud vertical extent h is possible. A second estimate of h comes from the correlation between the angular standard deviation of POLDER oxygen pressure σPO2 and the cloud vertical extent. This correlation is studied in detail in the case of liquid clouds. It is shown to be spatially and temporally robust, except for clouds above land during winter months. The analysis of the correlation's dependence on the scene's characteristics leads to a parameterization providing h from σPO2. For liquid water clouds above ocean in 2008, the mean difference between the actual cloud vertical extent and the one retrieved from σPO2 (from the pressure difference) is 5 m (−12 m). The standard deviation of the mean difference is close to 1000 m for the two methods. POLDER estimates of the cloud geometrical thickness obtain a global score of 50% confidence for a relative error of 20% (40%) of the estimate for ice (liquid) clouds over ocean. These results need to be validated outside of the CALIPSO/CloudSat track

Cloud thermodynamic phase inferred from merged POLDER and MODIS data

International audienceThe global spatial and diurnal distribution of cloud properties is a key issue for understanding the hydrological cycle, and critical for advancing efforts to improve numerical weather models and general circulation models. Satellite data provides the best way of gaining insight into global cloud properties. In particular, the determination of cloud thermodynamic phase is a critical first step in the process of inferring cloud optical and microphysical properties from satellite measurements. It is important that cloud phase be derived together with an estimate of the confidence of this determination, so that this information can be included with subsequent retrievals (optical thickness, effective particle radius, and ice/liquid water content). In this study, we combine three different and well documented approaches for inferring cloud phase into a single algorithm. The algorithm is applied to data obtained by the MODIS (MODerate resolution Imaging Spectroradiometer) and POLDER3 (Polarization and Directionality of the Earth Reflectance) instruments. It is shown that this synergistic algorithm can be used routinely to derive cloud phase along with an index that helps to discriminate ambiguous phase from confident phase cases. The resulting product provides a semi-continuous confidence index ranging from confident liquid to confident ice instead of the usual discrete classification of liquid phase, ice phase, mixed phase (potential combination of ice and liquid particles), or simply unknown phase clouds. This approach is expected to be useful for cloud assimilation and modeling efforts while providing more insight into the global cloud properties derived from satellite data

Cloud information content analysis of multi-angular measurements in the oxygen A-band: application to 3MI and MSPI

The vertical distribution of cloud cover has a significant impact on a large number of meteorological and climatic processes. Cloud top altitude and cloud geometrical thickness are then essential. Previous studies established the possibility of retrieving those parameters from multi-angular oxygen A-band measurements. Here we perform a study and comparison of the performances of future instruments. The 3MI (Multi-angle, Multi-channel and Multi-polarization Imager) instrument developed by EUMETSAT, which is an extension of the POLDER/PARASOL instrument, and MSPI (Multi-angles Spectro-Polarimetric Imager) develoloped by NASA's Jet Propulsion Laboratory will measure total and polarized light reflected by the Earth's atmosphere–surface system in several spectral bands (from UV to SWIR) and several viewing geometries. Those instruments should provide opportunities to observe the links between the cloud structures and the anisotropy of the reflected solar radiation into space. Specific algorithms will need be developed in order to take advantage of the new capabilities of this instrument. However, prior to this effort, we need to understand, through a theoretical Shannon information content analysis, the limits and advantages of these new instruments for retrieving liquid and ice cloud properties, and especially, in this study, the amount of information coming from the A-Band channel on the cloud top altitude (CTOP) and geometrical thickness (CGT). We compare the information content of 3MI A-Band in two configurations and that of MSPI. Quantitative information content estimates show that the retrieval of CTOP with a high accuracy is possible in almost all cases investigated. The retrieval of CGT seems less easy but possible for optically thick clouds above a black surface, at least when CGT > 1–2 km

An algorithm to retrieve ice water content profiles in cirrus clouds from the synergy of ground-based lidar and thermal infrared radiometer measurements

The algorithm presented in this paper was developed to retrieve ice water content (IWC) profiles in cirrus clouds. It is based on optimal estimation theory and combines ground-based visible lidar and thermal infrared (TIR) radiometer measurements in a common retrieval framework in order to retrieve profiles of IWC together with a correction factor for the backscatter intensity of cirrus cloud particles. As a first step, we introduce a method to retrieve extinction and IWC profiles in cirrus clouds from the lidar measurements alone and demonstrate the shortcomings of this approach due to the backscatter-to-extinction ambiguity. As a second step, we show that TIR radiances constrain the backscattering of the ice crystals at the visible lidar wavelength by constraining the ice water path (IWP) and hence the IWC, which is linked to the optical properties of the ice crystals via a realistic bulk ice microphysical model. The scattering phase function obtained from the microphysical model is flat around the backscatter direction (i.e., there is no backscatter peak). We show that using this flat backscattering phase function to define the backscatter-to-extinction ratio of the ice crystals in the retrievals with the lidar-only algorithm results in an overestimation of the IWC, which is inconsistent with the TIR radiometer measurements. Hence, a synergy algorithm was developed that combines the attenuated backscatter profiles measured by the lidar and the measurements of TIR radiances in a common optimal estimation framework to retrieve the IWC profile together with a correction factor for the phase function of the bulk ice crystals in the backscattering direction. We show that this approach yields consistent lidar and TIR results. The resulting lidar ratios for cirrus clouds are found to be consistent with previous independent studies.</p

Liquid cloud optical property retrieval and associated uncertainties using multi-angular and bispectral measurements of the airborne radiometer OSIRIS

In remote sensing applications, clouds are generally characterized by two properties: cloud optical thickness (COT) and effective radius of water–ice particles (Reff), as well as additionally by geometric properties when specific information is available. Most of the current operational passive remote sensing algorithms use a mono-angular bispectral method to retrieve COT and Reff. They are based on pre-computed lookup tables while assuming a homogeneous plane-parallel cloud layer. In this work, we use the formalism of the optimal estimation method, applied to airborne near-infrared high-resolution multi-angular measurements, to retrieve COT and Reff as well as the corresponding uncertainties related to the measurement errors, the non-retrieved parameters, and the cloud model assumptions. The measurements used were acquired by the airborne radiometer OSIRIS (Observing System Including PolaRization in the Solar Infrared Spectrum), developed by the Laboratoire d'Optique Atmosphérique. It provides multi-angular measurements at a resolution of tens of meters, which is very suitable for refining our knowledge of cloud properties and their high spatial variability. OSIRIS is based on the POLDER (POlarization and Directionality of the Earth's Reflectances) concept as a prototype of the future 3MI (Multi-viewing Multi-channel Multi-polarization Imager) planned to be launched on the EUMETSAT-ESA MetOp-SG platform in 2024. The approach used allows the exploitation of all the angular information available for each pixel to overcome the radiance angular effects. More consistent cloud properties with lower uncertainty compared to operational mono-directional retrieval methods (traditional bispectral method) are then obtained. The framework of the optimal estimation method also provides the possibility to estimate uncertainties of different sources. Three types of errors were evaluated: (1) errors related to measurement uncertainties, which reach 6 % and 12 % for COT and Reff, respectively, (2) errors related to an incorrect estimation of the ancillary data that remain below 0.5 %, and (3) errors related to the simplified cloud physical model assuming independent pixel approximation. We show that not considering the in-cloud heterogeneous vertical profiles and the 3D radiative transfer effects leads to an average uncertainty of 5 % and 4 % for COT and 13 % and 9 % for Reff.</p

POLDER observations of cloud bidirectional reflectances compared to a plane-parallel model using the International Satellite Cloud Climatology Project cloud phase functions

International audienceThis study investigates the validity of the plane-parallel cloud model and in addition the suitability of water droplet and ice polycrystal phase functions for stratocumulus and cirrus clouds, respectively. To do that, we take advantage of the multidirectional viewing capability of the Polarization and Directionality of the Earth's Reflectances (POLDER) instrument which allows us to characterize the anisotropy of the reflected radiation field. We focus on the analysis of airborne-POLDER data acquired over stratocumulus and cirrus clouds during two selected flights (on April 17 and April 18, 1994) of the European Cloud and Radiation Experiment (EUCREX'94) campaign. The bidirectional reflectances measured in the 0.86 μm channel are compared to plane-parallel cloud simulations computed with the microphysical models used by the International Satellite Cloud Climatology Project (ISCCP). Although clouds are not homogeneous plane-parallel layers, the extended cloud layers under study appear to act, on average, as a homogeneous plane-parallel layer. The standard water droplet model (with an effective radius of 10 μm) used in the ISCCP analysis seems to be suitable for stratocumulus clouds. The relative root-mean-square difference between the observed bidirectional reflectances and the model is only 2%. For cirrus clouds, the water droplet cloud model is definitely inadequate since the rms difference rises to 9%; when the ice polycrystal model chosen for the reanalysis of ISCCP data is used instead, the rms difference is reduced to 3%

Cloud thermodynamic phase inferred from merged POLDER and MODIS data

The global spatial and diurnal distribution of cloud properties is a key issue for understanding the hydrological cycle, and critical for advancing efforts to improve numerical weather models and general circulation models. Satellite data provides the best way of gaining insight into global cloud properties. In particular, the determination of cloud thermodynamic phase is a critical first step in the process of inferring cloud optical and microphysical properties from satellite measurements. It is important that cloud phase be derived together with an estimate of the confidence of this determination, so that this information can be included with subsequent retrievals (optical thickness, effective particle radius, and ice/liquid water content). In this study, we combine three different and well documented approaches for inferring cloud phase into a single algorithm. The algorithm is applied to data obtained by the MODIS (MODerate resolution Imaging Spectroradiometer) and POLDER3 (Polarization and Directionality of the Earth Reflectance) instruments. It is shown that this synergistic algorithm can be used routinely to derive cloud phase along with an index that helps to discriminate ambiguous phase from confident phase cases. The resulting product provides a semi-continuous index ranging from confident liquid to confident ice instead of the usual discrete classification of liquid phase, ice phase, mixed phase (potential combination of ice and liquid particles), or simply unknown phase clouds. The index value provides simultaneously information on the phase and the associated confidence. This approach is expected to be useful for cloud assimilation and modeling efforts while providing more insight into the global cloud properties derived from satellite data

Ice crystal shapes in cirrus clouds derived from POLDER-1/ADEOS-1.

International audienceThis paper discusses the retrieval of ice crystal shapes of cirrus clouds on a global scale using observations collected with POLDER-1 (POLarization and Directionality of the Earth Reflectance) onboard the ADEOS-1 platform. The retrieval is based on polarized bidirectional observations made by POLDER. First, normalized polarized radiances are simulated for cirrus clouds composed of ice crystals that differ in shape and are randomly oriented in space. Different values of cloud optical depths, viewing geometries and solar zenith angles are used in the simulations. This sensitivity study shows that the normalized polarized radiance is highly sensitive to the shape of the scatterers for specific viewing geometries, and that it saturates after a few scattering events, which makes it rapidly independent of the optical depth of the cirrus clouds. Next, normalized polarized radiance observations obtained by POLDER have been selected, based on suitable viewing geometries and on the occurrence of thick cirrus clouds composed of particles randomly oriented in space. For various ice crystal shapes these observations are compared with calculated values pertaining to the same geometry, in order to determine the shape that best reproduces the measurements. The method is tested fully for the POLDER data collected on January 12, 1997. Thereafter, it is applied to six periods of 6 days of observations obtained in January, February, March, April, May, and June 1997. This study shows that the particle shape is highly variable with location and season, and that polycrystals and hexagonal columns are dominant at low latitudes, whereas hexagonal plates occur more frequently at high latitudes

Overview of the Chemistry-Aerosol Mediterranean Experiment/Aerosol Direct Radiative Forcing on the Mediterranean Climate (ChArMEx/ADRIMED) summer 2013 campaign

The Chemistry-Aerosol Mediterranean Experiment (ChArMEx; http://charmex.lsce.ipsl.fr) is a collaborative research program federating international activities to investigate Mediterranean regional chemistry-climate interactions. A special observing period (SOP-1a) including intensive airborne measurements was performed in the framework of the Aerosol Direct Radiative Impact on the regional climate in the MEDiterranean region (ADRIMED) project during the Mediterranean dry season over the western and central Mediterranean basins, with a focus on aerosol-radiation measurements and their modeling. The SOP-1a took place from 11 June to 5 July 2013. Airborne measurements were made by both the ATR-42 and F-20 French research aircraft operated from Sardinia (Italy) and instrumented for in situ and remote-sensing measurements, respectively, and by sounding and drifting balloons, launched in Minorca. The experimental setup also involved several ground-based measurement sites on islands including two ground-based reference stations in Corsica and Lampedusa and secondary monitoring sites in Minorca and Sicily. Additional measurements including lidar profiling were also performed on alert during aircraft operations at EARLINET/ACTRIS stations at Granada and Barcelona in Spain, and in southern Italy. Remote-sensing aerosol products from satellites (MSG/SEVIRI, MODIS) and from the AERONET/PHOTONS network were also used. Dedicated meso-scale and regional modeling experiments were performed in relation to this observational effort. We provide here an overview of the different surface and aircraft observations deployed during the ChArMEx/ADRIMED period and of associated modeling studies together with an analysis of the synoptic conditions that determined the aerosol emission and transport. Meteorological conditions observed during this campaign (moderate temperatures and southern flows) were not favorable to producing high levels of atmospheric pollutants or intense biomass burning events in the region. However, numerous mineral dust plumes were observed during the campaign, with the main sources located in Morocco, Algeria and Tunisia, leading to aerosol optical depth (AOD) values ranging between 0.2 and 0.6 (at 440 nm) over the western and central Mediterranean basins. One important point of this experiment concerns the direct observations of aerosol extinction onboard the ATR-42, using the CAPS system, showing local maxima reaching up to 150Mm(-1) within the dust plume. Non-negligible aerosol extinction (about 50Mm(-1)) has also been observed within the marine boundary layer (MBL). By combining the ATR- 42 extinction coefficient observations with absorption and scattering measurements, we performed a complete optical closure revealing excellent agreement with estimated optical properties. This additional information on extinction properties has allowed calculation of the dust single scattering albedo (SSA) with a high level of confidence over the western Mediterranean. Our results show a moderate variability from 0.90 to 1.00 (at 530 nm) for all flights studied compared to that reported in the literature on this optical parameter. Our results underline also a relatively low difference in SSA with values derived near dust sources. In parallel, active remote-sensing observations from the surface and onboard the F-20 aircraft suggest a complex vertical structure of particles and distinct aerosol layers with sea spray and pollution located within the MBL, and mineral dust and/or aged North American smoke particles located above (up to 6–7 km in altitude). Aircraft and balloon-borne observations allow one to investigate the vertical structure of the aerosol size distribution showing particles characterized by a large size (> 10 μm in diameter) within dust plumes. In most of cases, a coarse mode characterized by an effective diameter ranging between 5 and 10 μm, has been detected above the MBL. In terms of shortwave (SW) direct forcing, in situ surface and aircraft observations have been merged and used as inputs in 1-D radiative transfer codes for calculating the aerosol direct radiative forcing (DRF). Results show significant surface SW instantaneous forcing (up to (-90)Wm(-2) at noon). Aircraft observations provide also original estimates of the vertical structure of SW and LW radiative heating revealing significant instantaneous values of about 5 K per day in the solar spectrum (for a solar angle of 30 ) within the dust layer. Associated 3-D modeling studies from regional climate (RCM) and chemistry transport (CTM) models indicate a relatively good agreement for simulated AOD compared with observations from the AERONET/PHOTONS network and satellite data, especially for long-range dust transport. Calculations of the 3-D SW (clear-sky) surface DRF indicate an average of about -10 to -20Wm(-2) (for the whole period) over the Mediterranean Sea together with maxima (-50Wm(-2)) over northern Africa. The top of the atmosphere (TOA) DRF is shown to be highly variable within the domain, due to moderate absorbing properties of dust and changes in the surface albedo. Indeed, 3-D simulations indicate negative forcing over the Mediterranean Sea and Europe and positive forcing over northern Africa. Finally, a multiyear simulation, performed for the 2003 to 2009 period and including an ocean–atmosphere (O–A) coupling, underlines the impact of the aerosol direct radiative forcing on the sea surface temperature, O–A fluxes and the hydrological cycle over the Mediterranean.French National Research Agency (ANR) ANR-11-BS56-0006ADEMEFrench Atomic Energy CommissionCNRS-INSU and Meteo-France through the multidisciplinary programme MISTRALS (Mediterranean Integrated Studies aT Regional And Local Scales)CORSiCA project - Collectivite Territoriale de Corse through Fonds Europeen de Developpement Regional of the European Operational ProgramContrat de Plan Etat-RegionEuropean Union's Horizon 2020 research and innovation program 654169Spanish Ministry of Economy and Competitivity TEC2012-34575Science and Innovation UNPC10-4E-442European Union (EU)Department of Economy and Knowledge of the Catalan Autonomous Government SGR 583Andalusian Regional Government P12-RNM-2409Spanish Government CGL2013-45410-R 26225